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ABSTRACT 

We present the first high-angular resolution (up to 0'/7, ^5000 AU) polarization and thermal dust 
continuum images toward the massive star-forming region W51 North. The observations were carried 
out with the Submillimeter Array (SMA) in both the subcompact (SMA-SubC) and extended (SMA- 
Ext) configurations at a wavelength of 870 //m. W51 North is resolved into four cores (SMA1 to 
SMA4) in the 870 /im continuum image. The associated dust polarization exhibits more complex 
structures than seen at lower angular resolution. We analyze the inferred morphologies of the plane- 
of-sky magnetic field (Bj_) in the SMA1 to SMA4 cores and in the envelope using the SMA- Ext and 
SMA-SubC data. These results are compared with the B^ archive images obtained from the CSO 
and JCMT. The polarization percentage is about 1 % to 4 %, and it is found to decrease with higher 
intensity in our SMA images, which is a similar trend as previously reported in the CSO and JCMT 
data. A correlation between dust intensity gradient position angles (</>vi) and magnetic field position 
angles (<^>b) is found in the CSO, JCMT and both SMA data sets. This correlation is further analyzed 
quantitatively. A systematically tighter correlation between </>vi and 4>b is found in the cores, whereas 
the correlation decreases in outside-core regions. Magnetic field-to-gravity force ratio (Eb) maps ar e 
derived using the newly developed polarization - intensity gradient method by iKoch et all (|2012aft . 
We find that the force ratios tend to be small (Sb 0.5) in the cores in all 4 data sets. In regions 
outside of the cores, the ratios increase or the field is even dominating gravity (Eb > 1)- This possibly 
provides a physical explanation of the tightening correlation between <^vi & n d <^b in the cores: the 
more the B field lines are dragged and aligned by gravity, the tighter the correlation is. Finally, we 
propose a schematic scenario for the magnetic field in W51 North to interpret the four polarization 
observations at different physical scales. 

Subject headings: ISM: individual (W51 d) - individual (W51 North) - individual (W51 IRS2) - 
individual (W51 A) - ISM: magnetic fields - polarization - stars: formation 



1. INTRODUCTION 

Stars form in giant molecul ar clouds under the threads 
of turbulence (see review by IMac Low fc Kless"en 12004ft 
and large-scale magnetic (B) fields (e.g. lCrutcherI [l999). 
Theoretically, the significance of the B field influences 
how structures ar e formed, such as the density contrast 
withi n structures (iKowal et al. 2007), the star formation 
rate (|Vazquez-Semadeni et al.1 12005: Price & Batel l2008r) 
and the suppressed fra gmentation ([Price fc Bate! I20071 
iHennebelle et all 120111 ). However, the B fields in star- 
forming clouds are not well constrained observationally, 
because they are difficult to detect. The Zeeman effect on 
spectral lines can trace the B field strength along the line 
of sight (e.g. iFjsh fc Reidl [20061: iVlemmings et all [20061: 
iFalgarone et al.l 120081: iSurcis et al.l 120121) . The B field 
properties on the plane of sky, B^, rely on polarization 
studies. 

Dust grains are thought to be partly aligned with 
their minor axes to the magnetic field lines. Several 
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alignment mec hanisms have been studied (see review by 
ILazar ian 2000) , and radiative tor ques are likely to be re- 



spons i ble fo r the dust alignme nt (Drainc & Wcingartncr 
19961 119971 : iCho fc Lazarianl 120051: ILazarian fc Hoan~ 
20071 ). Due to the difference in emissivity along the 



major and minor axes of the elongated dust grains, 
the thermal emission from dust grains can therefore 
be pa rtly linearly polar ized perpendicular to the field 
lines (Hildcbrand 1981). The Bj^ morphology can 
therefore be inferred by rotating the detected po- 
larization of the thermal emission by 90°. Linear 
polarization of thermal dust emission has been im- 
aged at far infrared and at submillim eter (submm) 

in s everal star-forming clouds fe.g.. | Greaves et al.l 

1994tlChrvsostomou et al.ll2002tlVaillancourt et al.ll2008r 
Dotson et all 120101: iMatthews et al.l I2009D . To further 



resolve the star-forming cores, higher angular resolution 
observations of the millimeter (mm) and submm con- 
tinuum emission are powerful, because the continuum is 
mostly optically thin at these wavelengths. Both mas- 
sive and low mass star-forming cores have been imaged 
in linear polarization with an angular resolutio n of a few 
arcsecond in the mm and submm regime (e .g . . lLai et al 
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;ond m trie mm and su bmm regime (e.g. . lLai et al 
iGirart et al.1 120061 [Tang et all l2009aH bl: iRao et al 



2009l : lTang et al.ll2010f ). W51 North is another good can 



didate to study the massive star formation processes. It 
has strong dust emission at mm wavelengths. This in- 
creases the possibilities to detect polarization, because 
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the linear polarization percentage is typically only a few 
percent in star-forming regions and therefore, difficult to 
measure. 

The W51 North region is active in star formation. It 
lies within the HII region complex G49.5-0.4 in W51 A, 
west of the W51 giant molecular cloud. It is located 
5— 8 kpc away in the Sagittarius spiral arm (jGenzel et al.l 
H98lllmai et all[200l Ku et al.ll2009fl . In this paper, we 
adopt a distance of 7 kpc. At radio wavelengths (1.3 
cm and 3.6 cm), an edge-brightened cometary HII re- 
gion, W5 1 d, and an ultraco mpact HII (UCHII) region, 
W51 d2, (jGaume et al.lll993t ) were detected The bright- 
est mm source is 2" to the east of W51 d2 (|Zhang et al.l 
U998tlSollins et al.ll2004HZapata et al.ll2008ft . Associated 
with this mm source, a group of H2O, OH and SiO masers 
were detected, and it has been also named as "the dom- 
inant center" (Schncps et al.l Il981[ ). In t his dominant 
center, infall signatures have be en detected (jSollins et al.1 
[200l iZaoata et all [20081 120091 ) . The outflows from the 
dominant center were r eported to have position angles 
between 105° and 150° (jEisner et al.ll2002t iZhang et all 
1998; Zapata et aill2009f) . An emb edded protostar candi- 
date OKYM 1 (also called KJD 3. lKraemer et al.ll200lD . 
1" to the north-east of the dominant center , was reported 
(jOkamoto et al.ll200i1: IFarbosa et al.ll200l . This region 
is, thus, a very active star-forming complex where several 
regions are in different evolutionary stages. 

Linear polarization of thermal dust continuum toward 
W51N was detected in the submm regime with single 
dish telescopes. At wav elengths of 35 /im with an an- 
gular resolution of 20", iDotson et al.l (|2010D show that 
the polarization is mainly u niform across W51 No r th. In 
cont rast, at 850 p,m both iChrvsostomou et"al"1 (|2002f ) 
and IMatthews et a l. (2009) show that there are more 
variations in the polarization orientations. Since there 
is probably multiple ongoing star formation in W51 
North, higher angular resolution measurements of the 
polarization are important to reveal the physical condi- 
tions at this scale. Her eafter, the data obtained f rom 
IDotson et all {2010) and IChrvsostomou et all (|2002D are 
called CSO and JCMT, respectively. 

Here, we report the new Submillimcter Array (SMA) 
870 /im continuum images and the detection of linearly 
polarized emission of dust continuum at an angular res- 
olution up to 0'.'7 toward the W51 North. Polarization 
percentages and position angles, and the correlation be- 
tween position angle and intensity are studied. With the 
goal of providing a more complete scenario of the role of 
the magnetic field over several scales, we further compare 
our new polarization results at a scale of 60 mpc with the 
published polarization images at a scale of a few pc at 
the submm regime from CSO and JCMT. The observa- 
tions and data reduction are described in section [2] The 
observational results are shown in section [3] In section @J 
we will present the detailed quantitative analyses of the 
polarization of our new observational results, and we will 
provide a physical explanation of the derived values. We 
then propose a schematic scenario for the magnetic field 
in W51 North. We draw conclusions and summarize in 
section [SJ 



2. OBSERVATIONS AND DATA REDUCTION 



TABLE 1 

Observational parameters 



Parameter 


SMA-Ext 


SMA-SubC 


Date 


2008 Jul 13 


2009 Sep 11 
6 


Available Antennae 


7 


Gain Calibrator 


1751+096, 1925+211 


1751+096 


Flux Calibrator 


Titan 


Uranus 


Bandpass/Pol Calib. 


3c454.3 


3c454.3 


Baseline range (kX) 


30 to 262 


10 to 30 


<Jl (Jy beam -1 ) 


0.06 


0.19 


<Ti (mjy beam -1 ) 


3 


15 


^syn 


Of.' 7 X Of.' 6 


4('lx3'/6 


PA. of 6» syn 


-57° 


20° 



Note. — a\ and o\ are the noise levels of Stokes I and 



the polarized intensity (ip), respectively. 

The observat ions were carried out using the SMA 
(|Ho et al.ll2004f) 5 in the extended and subcompact con- 
figuration, respectively. In the following, the results ob- 
tained from the extended and the subcompact configu- 
ration are called SMA-Ext and SMA-SubC, respectively. 
The synthesized beams (angular resolutions), 9 syn , of 
SMA-Ext and SMA-SubC are - 0"7 and 4", respectively. 
In both observations, the local oscillator frequency was 
tuned to 341.482 GHz (870 /mi). The observational de- 
tails are listed in Table 1. 

The SMA-Ext observation was carried out along with 
W51 e2/e8. These two sources shared the same cali- 
brators and were reduced and imaged in the same way. 
The details of the observation and data reduction are 
described in Tang et al. (2009). In brief, the SMA- 
Ext data were reduced and imaged in MIRIAD following 
standard procedures. An additional calibration on the 
instrumental polarization (also called the leakage terms) 
was derived using the calibrator 3c454.3. The largest size 
scale which could be sampled in the extended track was 
~8" (0.3 pc). For the subcompact track, the data were 
reduced with MIR and then imaged using MIRIAD. 

The position angle of the polarization, <p po i, is cal- 
culated from: <f> po i = \ tan -1 ^, increasing counter- 
clockwise in the range —90° to 90°. The strength, I p , 
and percentage, P(%), of the linearly polarized emission 
are calculated from: I 2 = Q 2 + U 2 - Oq u and P{%) — 
I p /I, respectively. Here, <7q,u is the noise level of the 
Stokes Q and U images, and it is also the bias correction 
due to the positive measure of I p (Leahy 1989; Wardle 
& Kronberg 1974). 

The primary beam (field of view) of the SMA at 345 
GHz is —30". All of the presented SMA maps in the fol- 
lowing have been corrected for the primary beam atten- 
uation. The presented maps were constructed with nat- 
ural weighting. The phase center is at Right Ascension 
(J2000)=19 /l 23 ro 40.05 s , Declination (J2000)=14°31'5 / /0. 

3. RESULTS 
3.1. Continuum Emission 
3.1.1. Dense Structures: SMA 1 to SMA4 

5 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



Dust continuum and polarization in W51 North 
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Fig. 1. — (a): C ontinuum map at 870 /im obtained from SMA-Ext and continuum map at 7 mm from EVLA (image adopted from 
IZapata et alj [20081) shown in white contours and c olor scale, respectiv ely. The small squares mark the position of the edge-brightened 
cometary /shell-like HII region W51 d and W51 d2 (Gaumc ct al. 1993|). The infrared source KJD3 (also called OKMY 1) is marked as *. 
The stars mark the identified O stars by Okamoto et al. (2001). The pluses mark the resolved submillimeter sources SMA 1 to 4. (b) and 
(c): Polarization (segments) and polarized intensity (grey scale) from the subcompact track (panel b) and extended track (panel c). The 
lengths of the polarization segments are proportional to the polarization percentages. The plotted thick and thin segments are above 3 and 
between 2 to 3 <ri p , respectively. The synthesized beams are shown as ellipses in the lower-left corners. Contours are at -3 (dash) and 3, 6, 

10, 20, 30, 40, ... X 0.06 and 0.19 Jy beam" 1 for SMA-Ext and SMA-SubC, respectively. 



TABLE 2 

Dust components 



Name I R.A. Dec. Flux (Jy) Deconvolved Size (" X "), PA(°) (<t> po \) Note 

SMA1 I 19:23:40.047 14:31:5.53 8.69±0.27 (0.98±0.03) X (0.68±0.02), -89±3 ^25 51 dominant center 

SMA2 19:23:39.970 14:31:5.70 2.59±0.24 (1.12±0.07) X (0.28±0.11), 39±3 4 30 

SMA3 19:23:39.828 14:31:5.14 2.41±0.27 (1.07±0.09) X (0.59±0.09), -88±7 -11 61 W51d2 

SMA4 19:23:39.756 14:31:5.37 1.40±0.22 (0.79±0.12) X (0.49±0.13), -23±21 -30 24 

Note. — Resolved components of the east- west dust ridge in the SMA-Ext data. (</> po i) and cr^ o] are the mean and 
the standard deviation of the polarization position angles (</> po i), respectively. Deconvolved sizes for SMA1 to SMA4 are 
obtained by fitting 2D-Gaussians. 



The 870 ^tm continuum emission is resolved in both 
SMA-SubC and SMA-Ext. With 6» syn of 4" (SMA-SubC), 
the 870 /im continuum emission is elongated in the east- 
west (EW) direction with an extension toward the north 
(Figurcf!} )). This is cons istent with the 2 mm continuum 
map (Zha ng et aLlll998l) and the thermal NH 3 emission 
(jZhang fc Hoi 119951 ). With 9 syn of 0'.'7 (SMA-Ext), the 
EW ridge is further resolved into SMA1 to SMA4, as la- 
beled in Figure [T^,. The coordinates and flux densities of 
SMAl to SMA4 are determined by fitting 2-D Gaussians. 
They are listed in Tabic H 

The brightest 870 /mi continuum peak is SMAl, which 
is consi stent with the loca tion of the brig htest peak at 
7 mm ([Zapata et al.l [20091 ) and at 2 mm (jZhang et al.l 
[19981 ). SMA2, SMA3 and SMA4 are resolved for the 
first time in the mm/submm regime. For SMA2, al- 
though there is no previous re port of detection, i t can 
be seen at 1.3 mm (Figure 1 in lZapata et~aT] |2008). As- 
sociated with SMA3 is the UC HII region W51 d2, 0"5 to 
the south, and a NBU maser (iMauersberger et al.llT987t 
iWilson et all 119911: iGaume et al.lll993l) . Among these 4 
resolved SMA sources, SMA3 is the most evolved source 
due to its free-free emission, which suggests the existence 
of ionized gas. The existence of masers in SMAl and 
SMA3 suggests that these two locations are active in star 
formation. 



3.1.2. Flux Density and Mass 

The flux density detected with SMA-SubC is 27±3 Jy. 
Comparing to the single dish measurement at the same 
location of 44 Jy (Chrysostomou et al. 2002), about 60% 
of the flux density is recovered in SMA-SubC. The com- 
pact structures (SMAl to SMA4) identified from SMA- 
Ext have a total flux density of 15±1 Jy over the EW 
dust ridge. Comparing to the total emission of 44 Jy, 
these compact structures make for about 34% of the to- 
tal emission. The majority of the emission is, thus, in the 
extended structure which is filtered out by the sampling 
of longer baselines. 

The free-free continuum emission is 2.5 mJy at 3.6 
cm toward W51 d 2 (named as SMA3 in this paper; 
IGaume et al.| [l993). and it will be 1.7 mJy or 16 mJy 
when extrapolated to 870 /im, assuming an optically thin 
spectral index (a) of —0.1 or an optically thick a of 0.6, 
respectively. While the flux density of SMA3 is 2.4 Jy at 
870/mi, the contribution from free- free emission is neg- 
ligible. The 870^m continuum emission along SMAl to 
SMA4 is, thus, mainly from thermal dust emission. For 
the extended HII region W51 d, the flux density is 1 Jy 
at 1.3 cm, and it will be 0.8 Jy or 5 Jy at 870 /im assum- 
ing a of —0.1 or 0.6, respectively. At this location, the 
detected intensity at 870 /im with SMA-SubC is ^5 Jy 
beam -1 . Thus, the traced extension to the north at 870 
/im could have a significant contribution from the free- 
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free emission, if this free-free emission is not optically 
thin between the wavelength of 1.3 cm and 870 fim. 

The total gas mass, M gas , can be estimated from 
the du st continuum, ass uming a gas temperature of 
100 K (jZhang et al.lll998T ). a dust opacity «A=87u ( um <x 
A- 2 « 0.4 cm g (following the standard equation 
in lAgladze et al.l [19 96) . and a gas-to-dust ratio of 100. 
For SMA1, M gas is ~130 roughly consistent with 
the value rep o rted i n lZhang et all (| 1998ft of 100 M Q and 
IZaoata et all IpOOfift of 90 M . For SMA 2, SMA 3 and 
SMA 4, M gas is 39, 36 and 21 M , respectively. We note 
that the estimated mass can vary by a factor of a few due 
to the uncertainty in k\, the assumed dust temperature, 
the gas-to-dust ratio and the distance to the source. 

3.2. Dust Polarization 
3.2.1. Morphology 

The linear polarization is detected in most parts of the 
dust ridge (Figures [TJd and Hh) . The presented polariza- 
tion segments were gridded every 0'.'3 for SMA-Ext and 
every 2" x l'/8 for SMA-SubC, which is about half of the 
synthesized beam in order to present the variations of 
the position angles (P.A.s). The detected polarization 
data above 2.5 <jj for the SMA-Ext and above 2 07 for 
the SMA-SubC are listed in Table H We choose to set 
a 2cr/ p threshold for the SMA-SubC in order to recover 
more polarization segments, where otherwise only 8 (out 
of 15) for SMA-SubC and 43 (out of 48) for SMA-Ext 
segments would remain with a 3a 1 limit. 

In the SMA-SubC map (Figure [TJ)) , there is a region 
without significant polarization about 2" north of SMA3. 
As zooming in with SMA-Ext, there are zones without 
significant polarization in between SMA1 and SMA2, be- 
tween SMA2 and SMA3 and at the peak of SMA3 and 
SMA4. The areas without significant polarization near 
the intensity peaks and in between SMA1 and SMA2 may 
be due to complex underlying field geometries, as they 
appear in between two dramatically different PA regions. 
One example of a complex field in a region with low/no 
polarization detection is the W51 e2/e8 region. Obser- 
vationally, the W51 e2 core was found to have low/no 
polarization in the 1.3 mm continuum with lower angu- 
lar resolution (3") bv lLai et al.l (|2001ft . This core has 
actually been detected and resolved in pol arization at 
870 11m at higher angular resolution (0.7") bv lTang et al.l 
(|2009bft . 

3.2.2. Distribution Plots 

The distributions of the polarization position angles 
(<^poi) are shown in Figure [2ji. In the shorter baseline 
data (SMA-SubC), the 4> vo \ are mainly in two ranges, 
namely around 50° and around -60°. As we will discuss 
in section 14. 1[ this bimodal-like distribution reflects a 
mirror symmetry in the polarization structure. The com- 
ponent at 50° is consistent with the single dish polariza- 
tion measurement at 350 /im with an angular resolution 
of about 20" (Dotson et al. 2010). With smaller 6> syn 
(SMA-Ext), the dust ridge and polarization are further 
resolved and exhibit more structures. In the majority 
of the detected patches, values for <f> vo \ are around —50° 
with smooth variations within these areas (Figure [TJ;). 
The distribution for SMA-Ext is, consequently, peaking 
around —50° in Figure^. Near SMA2, the po i are very 



TABLE 3 

Polarization data detected with the SMA 



AR.A. 


ADcc. 


P 


</>pol 


°>pol 


h 


1 


(") 


(") 


(%) 


(°) 


(°) 


(mjy/b) 


(Jy/b) 


Extended Track 












0.3 


-0.3 


1.8 ± 0.6 


-58 


10 


9 


0.48 


0.0 


-0.3 


2.2 ± 0.5 


-55 


6 


11 


0.62 


-0.3 


-0.3 


2.0 ± 0.6 


-52 


9 


9 


0.44 


-2.7 


-0.3 


3.2 ± 0.8 


-66 


7 


11 


0.36 


-3.0 


-0.3 


3.4 i 0.6 


-67 


5 


15 


0.43 


-3.3 


-0.3 


2.8 i 0.7 


-58 


7 


11 


0.39 


-3.6 


-0.3 


3.9 i 0.8 


-49 


6 


14 


0.36 


-3.9 


-0.3 


4.3 ± 1.0 


-47 


6 


12 


0.29 


0.3 


0.0 


0.7 ± 0.2 


-70 


10 


8 


1.13 


0.0 


0.0 


0.8 ± 0.2 


-57 


7 


11 


1.50 


-0.3 


0.0 


0.9 ± 0.2 


-47 


8 


10 


1.17 


-2.7 


0.0 


2.3 ± 0.6 


-97 


7 


12 


0.50 


-3.0 


0.0 


1.2 i 0.4 


-83 


8 


10 


0.80 


-3.6 


0.0 


1.4 i 0.4 


-19 


9 


9 


0.63 


-3.9 


0.0 
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- 11 


8 


11 


0.45 


-4.2 
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11 
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-81 


9 
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Note. — Offset in Right Ascension (AR.A.) and offset in Dec- 
lination (ADec.) with respect to the phase center (19 ,1 23 m 40.05 s , 
14°31'5'.'0)j2000 1 polarization percentage (P), position angle of po- 
larization ((fipoi), the uncertainty of the position angle of polariza- 
tion (f<^ pol ), polarization intensity (/ p ), and Stokes I intensity (/). 

I p and I are in units of mjy beam -1 (mjy/b) and Jy beam -1 
(Jy/b), respectively. Data listed are above 2.5<tt p and above 2(Tj p 
for SMA-Ext and SMA-SubC, respectively. 
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Fig. 2. — (a) Distributions of detected polarization position angles (</> po i). (b) Polarization percentage (/'%) versus normalized Stokes I 
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different from the po i at the other positions. The po i 
of these segments are more perpendicular to the inten- 
sity gradients at these positions (Figure [TJ;) . Since there 
are no reported ongoing star formation activities here, 
the origin of this abrupt change in po i from SMA1 to 
SMA2 and to SMA3 is not clear. The component with 
4> po \ of 50° seen in the SMA-SubC map is not seen in the 
SMA-Ext map (Figure [Tb). This suggests that this com- 
ponent originates from extended structures and is not 
sampled in the SMA-Ext data. 

A log-log plot of the polarization percentage (P) ver- 
sus the normalized intensity is shown in Figure^). P(%) 
is in the range of <1 to 4.5% and the mean P(%) is 
1.81% and 1.07% for the SMA-Ext and SMA-SubC, re- 
spectively. There is a general trend for both our high- 
angular-resolution data and the single dish data from 
the CSO and the JCMT: P decreases with larger inten- 
sity. Numerical simulations have shown that such anti- 
correlations can be due to a more complex DO i- structure 
in dense regions ([Falceta-Goncalves et al.ll2008t) . 

4. DISCUSSION 

In this section we discuss the different morphologies in 
the dust emission with its associated B field from 4 dif- 
ferent scales/layers. The larger scale data (shown in Fig- 
ures |3fe, and |3p) ar e from the si ngle dish by CSO ( Hertz, 
IDowell et all il998) at 350 //m (IDotson et all 120101) at a 
scale of 3 pc and by JCMT (SCUBAlHoUand et al.l ll998D 
at 850 /im (jChrvsostomou et alj|20020 at a scale of 2 pc. 
The smaller scale data are from our SMA-SubC at a scale 
of 0.3 pc and from SMA-Ext at a scale of 60 mpc at 870 
fjbrn (Figure [3b ) . The position angle of the magnetic field 
(0b) is obtained by rotating the detected <j) po \ by 90°. 
We also provide a quantitative analysis in support of the 
observational description in section |4~2"1 We conclude in 
section 14.31 with a schematic scenario of the magnetic 
field observations so far in the W51 North region. 

4.1. SMA high-resolution magnetic field morphologies 

The dust continuum morphology in the SMA-SubC 
shows regular contours with an extension to the west. 
These contours show some symmetries along an EW 



axis around ADec.=0" (Figure [3b) . Similarly, the de- 
tected polarization segments reveal a north-south sym- 
metry (upper and lower planes) along the same axis. As- 
suming flux-freezing with the gas mostly moving along 
the field lines, the combined observed symmetry pat- 
terns (in the dust continuum and the field morphologies) 
are suggestive to find mass accumulations mostly along 
the symmetry axis. Indeed, the detected cores (SMAl 
to SMA4) in the SMA-Ext data (Figure [3b) are mostly 
aligned on an EW axis at the above declination. Thus, 
the observed field configuration on the ~ 0.3 pc scale 
(SMA-SubC) appears to be channeling and aligning ma- 
terial, which is detected as denser cores on a ~ 60 mpc 
scale (SMA-Ext) once the extended emission is filtered 
out. 

Unlike the symmetric field configuration on the ~ 
0.3 pc scale, the smaller individual cores on the ~ 60 mpc 
scale show very different field morphologies. SMAl ex- 
hibits 4>b orientations roughly comparable to that of the 
one corresponding segment of the subcompact data (Fig- 
ure [5b). The morphology possibly also shows a hint of 
gravity dragging in the B field lines toward the center. 
The patch of segments in SMA2 is rather orthogonal to 
the SMA-SubC segment next to it. SMA3 and SMA4 
show again a more closer alignment with the single SMA- 
SubC segment but with some deviations up to about 45°. 
Thus, within one resolution element of the SMA-SubC 
data (where neighboring segments along the symmetry 
axis show only little changes in orientation) , the detected 
4>b in the SMA-Ext data show changes in orientation of 
up to 90°. 

The B morphology associated with SMAl is particu- 
larly interesting. As mentioned in section [TJ collapsing 
signatures and a rotati ng molecular ring have been re- 
ported tow ard SMAl (jSollins et al.l 120041: 1 Zapata et "all 
120081 12009D . The detected pinched B field morphology 
associated with SMAl is, therefore, consistent with the 
picture of the B field lines being dragged by inward mo- 
tions and/or the rotation of the molecular ring. 

Despite the very different field orientations in SMAl 
to SMA4, the </>b in these cores all reveal a common 
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Fig. 3. — The B field map in W51 North obtained from four different scales and/or layers (see also Figurc[7)l: 3 pc scale by HERTZ on 
the CSO in (a) and 2 pc scale by SCUBA on the JCMT in (b), and sub-pc scale by the SMA-SubC (black segments) and by SMA-Ext 
(white segments) in (c). In (a) and (b), the contours and the red segments label the continuum emission and the B field orientations, 
respectively. The small circles mark the data with larger uncertainties in position angles (15° to 30°) as compared to the other data points. 
The contours start from and step in 5% of the peak intensity detected in CSO of 694 Jy and in JCMT of 44 Jy in (a) and (b), respectively. 
The black square in panel (a) marks the mapping area in panel (b), and the one in panel (b) marks the mapping area in panel (c). In (c), 
the continuum emission traced with the SMA-SubC is shown in orange contours and in color scale, and the contours are 3, 6, 10, 20, 30 
X 0.19 (Jy beam -1 ). The continuum emission detected in the SMA-Ext is plot ted in grey contours in the same level as in Figure 1. The 
cyan X marks the center of the rotating toroid detected by Zapata ct al. (2009). 



trend: the field segments tend to be aligned with the 
orientations of the dust continuum intensity gradients, 
</>vi, ha many locations 6 . In section I4T21 we will have a 
more quantitative analysis on this correlation. 

4.2. Quantitative Analysis 

Figures [4^ to[4ji show the dust emission maps overlaid 
with the magnetic field segments (red) and the intensity 
gradient segments (blue) . Inspired by the close alignment 
of the position angles of the intensity gradients, </>vii and 
the position angles of the B field, 4>b, we apply further 
analyses in this subsection. 

4.2.1. Correlation of <f>\j\ and <f>B 

We first compare the correlation of (/>vi and </>b in these 
4 data sets. The correlation plots are shown in Figure 
The uniformity of the B field traced with the CSO is 
apparent here with a large group of data points scatter- 
ing around 4>b of 120°. Except for the CSO data, the 
correlation between 0vi and 4>b can be seen by eye. 

The correlation is further analyzed quantitatively using 
Pearson's correlation coefficient, C. The uncertainties in 

6 The dust intensity gradients are orthogonal to the dust emis- 
sion contours (in the plane of sky). The resulting position angles, 
</>Vl> are defined in the range of 0° to 180°. 



the correlation coefficient (the standard deviation), ac, 
are estimated using a Monte-Carlo approach where errors 
in both <^>b and </>vi are propagated through the calcu- 
lation. Gaussian distributions of the 0b (= po i + 90°) 
uncertainties at two sigma level, i.e. 2 times the <70 pol 
term in Table El are adopted. Gaussian distributions 
with mean uncertainties of ~ 3° — 5° - which is the re- 
sult after interpolating between Stokes / values in order 
to calculate gradients - are used for the c/!>vi uncertain- 
ties. Resulting values of ac are in the range of 0.02 to 
0.08. 

As listed in Table 21 C is highest in the highest resolu- 
tion data (SMA-Ext), being 0.88±0.01, suggesting that 
the correlation between 0yi and 4>b is tightest in the 
structures revealed with SMA-Ext. For SMA-SubC, the 
C value is 0.78±0.03. The C values for the CSO and 
the JCMT data are 0.81±0.03 and 0.70±0.02, respec- 
tively. We further derive the C values for the CSO and 
the JCMT data for two separated regions: the core and 
the region outside of the core. The C values are larger 
in the cores, being 0.88±0.04 and 0.73±0.02 in the CSO 
and JCMT data, respectively. This is similar to the high 
C value in SMA-Ext, further confirming that the </>vi 
and (f>B correlation is higher in the cores. Outside of the 
core region, the C values are generally smaller (0.59±0.04 
for CSO and 0.66±0.02 for JCMT), and the correlation 
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TABLE 4 

Analysis Summary 



Parameters 




CSO 






JCMT 




SMA-SubC 




SMA-Ext 




f(syn) (") 




20 






9.3 




4 




0.7 




A (^tm) 




350 






850 




870 




870 




B morphology 




uniform 




uniform patches, 


changes 


shaped, cometary-likc, 


shaped, possibly hourglass-like 










in larj 


;e scale 




N-S mirror symmetry 










all 


core" 


out" 


all 


core" 


out" 


all 


all 


corel c 


core2° 


C 


0.81 


0.88 


0.59 


0.70 


0.73 


0.66 


0.78 


0.88 


0.89 


0.86 




0.06 


0.08 


0.08 


0.03 


0.04 


0.03 


0.05 


0.02 


0.02 


0.03 


<|A<M) 


12° 






43° 






43° 


39° 






CT |A0 B l 


13° 






27° 






26° 


26° 






<l«l> 


40° 


51° 


34° 


44° 


38° 


48° 


37° 


35° 


34° 


36° 


y\ 


22° 


29° 


15° 


27° 


28° 


26° 


26° 


24° 


26° 


23° 


<£ B > 


0.71 


0.69(0.59) 


0.72 


1.17 


0.42(0.27) 


1.84 


0.47(0.18) 


0.48(0.38) 


0.47(0.37) 


0.48(0.39) 



Note. — Statistical quantities derived using the position angles of the intensity gradient (0vi) an d the position angles of the magnetic field 
(0b). #(syn) an d ^ are the angular resolution and the wavelength of the traced emission, respectively. C is the Pearson's correlation coefficient 
between 0b and 0vi, an d its standard deviation is given by ac ■ (|A0b|) and cT|a^ b I arc mean difference of 0b of polarization pairs and 
the standard deviation of the difference, respectively. and crj 5 is the mean of the absolute difference between 0b and 0vi an d its standard 

deviation, res pectively. (Sb) is th e mean of the ratio of the magnetic field force to gravitational force (Sb), which is calculated with the method 
introduced in Koch ct al. (2012a). Values after removing the outliers arc listed in parenthesis. 



( a ) Values derived using the pixels marked within the cyan square in Figure EH 
(k) Values derived using the pixels outside the cyan square in Figure ^ (named . 
( c ) Corel and core2 refers to SMA1,2 and SMA3,4, respectively. 



outside of the core region in the text). 



is worse. It is systematic that both the CSO and the 
JCMT data show a difference in C between the core and 
the outside-the-core region. We will propose a physical 
interpretation for such a change in C in section f4. 2.21 

In addition to the </>vi-</>b plot, we also compare the 
distribution of the angle difference, <fis, between </>b and 
0vi (Figure [6^). It is important to note that all the 
distributions are non-Gaussian, i.e., the differences in the 
magnetic field and intensity gradient orientations are not 



only due to a measurement uncertainty, but point toward 
a more fundamental p hysical interpretation as explored 
in IKoch et al.l (|2012ali . 

The differences of individual pairs of field segments 
over an observed map are also compared. The histograms 
in Figure[6jD display | Ac^b | = |</>B,i — 0b.j for all possible 
pairs of segments i and j in a map. A systematic differ- 
ence is apparent between the CSO data and the other 3 
data sets. The CSO distribution peaks at small values, 
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Fig. 5. — Position angles of the intensity gradient orientations, </>vii an d the B field orientations, ipB, in solid squares. Crosses mark the 
absolute differences between </>vi and <j)~B- In order to properly display their correlations, the 0b is re-defined beyond the 0°-180° range 
(filled squares above 180° and below 0°) for ambiguous situations; i.e., a pair of position angles of 175° and 5° is displayed as 175° and 
185°. The hallowed squares mark their original data points. Note that this re-definition does not change the correlation in any way, but is 
merely a consequence of the 180° ambiguity in the range of position angles. The straight line represents a perfect correlation. 



with an average change in orientation (|A0b|) = 12°, 
whereas the other distributions are close to fiat with av- 
erage changes around 40° (Table 2|) . This reflects the 
visual impression that the field is mostly uniform in ori- 
entation in the CSO map at of 120°. In the other 
maps, the field orientations seem to be more significantly 
shaped and overcome by other forces, erasing the original 
likely uniform orientation of the field. Thus, one might 
conclude that on the largest scales, a uniform field mor- 
phology points to a more significant role of the magnetic 
field. 

4.2.2. Application of the Polarization - Intensity Gradient 

Method 

In this subsection, we apply the new polarization- 
intensity gradient method to derive the force ratio (Eb) 
map of the magne tic field to the gravitational force (see 
IKoch et al.l I2012al for the development of the method). 
The derivation of Eb makes use of \4>s\ in combination 
with an additional angle between </>vi and the local grav- 
ity direction. This local gravity direction is calculated 
from the positions of the gravity centers which are iden- 
tified with the emission peaks in a map 7 . For SMA- 
Ext, two gravity centers are defined, namely SMA1 and 
SMA3, around which most of the polarization is detected. 
Figures to |3fr show maps of the force ratio E# ■ Val- 
ues averaged over the map, (Eb), are listed in Table |4j 
The uncertainties of the derived ratios depend on several 
factors, such as the 4>po\ uncertainties and the positions 
of the gravity centers in a map. Adopting uncertainties 
like in the previous section for C orr , typical errors in Eb 
are ~ 5% for core-like structures like in SMA-Ext. The 
errors can grow to ~ 20% for larger scale field morpholo- 
gies where gravity centers are less defined. Averaging 
Eb over a larger area (cores, outside-core regions) will 

7 For the origi nal d erivati on of the force ratio Eb we refer the 
reader to IKoch et al.l i2012aT ) . Eb is found from an MHD force 
equation, comparing the field tension force -Fb and the gravita- 
tional force Fq via measurable angles: Eb = Fb/Fq = sin ip/ sin a. 
The angle tp is defined by the intensity gradient and the direction of 
gravitational force. The angle a is between the intensity gradient 
and the originally detected polarization orientations. 



typically lead to errors of a few percents. 

The (Eb) values of the large-scale observations with 
the CSO and the JCMT are around 0.71 and 1.17 when 
averaged over the entire area, respectively. If we separate 
the values into core region and outside-the-core region, 
(Eb) is 0.69 (~ 0.59 if outliers are removed) and 0.72 for 
CSO in the core and outside, respectively. For JCMT, 
the variation is more dramatic, being 0.42 (0.27 if outliers 
are removed) and 1.84 for the core and the region outside 
of the core, respectively. We, thus, find that Eb tends to 
be smaller for the cores than in the areas outside of the 
cores. The smaller-scale SMA data also show a reduced 
magnetic field significance with (Eb) ~ 0.5. Generally, 
closer to and around the emission peaks, the force ratios 
seem to be smaller. This implies that on these scales the 
magnetic field is more shaped and dominated by gravity. 
In the regions outside of the core, the ratios are larger or 
the field is even dominating gravity (Eb > 1 for JCMT). 

This change in the role of the magnetic field based on 
the force ratio Eb might provide a physical explanation 
of the change of the correlation coefficient C described 
in section 14.2.11 The tighter correlation of </>vi and 0b 
in the cores suggests that the B field lines are mainly 
dragged by gravity and, consequently, align with the in- 
tensity gradient. Therefore, quantitatively, the B field 
is indeed of minor importance in the cores. We empha- 
size that the C value is purely statistical without any 
physical assumptions. In contrast, the Eb value is calcu- 
lated by identifying various terms in an ideal MHD force 
equation in an observed Stokes / and dust polarization 
map. These two approaches are, thus, largely indepen- 
dent. They, nevertheless, seem to lead to a consistent 
picture where the magnetic field significance is reduced 
in the core regions. 

4.3. Global Picture of the Magnetic Field in W51 North 

Figure[7]schematically summarizes the 4 magnetic field 
observations analyzed in the previous sections. Each 
panel represents a plane-of-sky projected observation of 
the CSO, JCMT, SMA-SubC and SMA-Ext, respec- 
tively, where a subsequently higher-angular-resolution 
observation is zooming in on a previous one. In this 
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Fig. 6. — (a) Histograms of the differences in position angles between the B field and the intensity gradient (tf>g). (b): Histograms of the 
absolute values of the differences in position angles of the polarization pairs (|A0b|)- 




sequence of observations, field morphologies are sampled 
from a 3 pc scale down to a 60 mpc physical scale. 

These four panels represent the field geometry from dif- 
ferent layers or from different scales in the W51 North re- 
gion. For the CSO image, the flux from an optically thick 
20" (the CSO resolution) source at a rotational temper- 
ature T of 15 K will be 270 Jy, and 760 Jy at T=30 K. 
This is in the range of Sto kes I intensity in the observed 
region (jDotson et al.ll2010h . We note that the rotational 
temperature of W5 1 N is 45 K as traced by ammonia lines 
(| Zhang fc Holll997f l and the detected structure is -10". 
A decreasing temperature from the UCHII region to the 
outer surrounding gas is further proposed bv lZhang et al.l 
(1998). Thus, the gas temperature within a 20" beam is 
unlikely to be much higher than 30 K. Unless the temper- 



atures are substantially higher, the average opacity must 
be of order 1. Furthermore, the observed polarization 
arises from the difference in emissivity of dust grains. In 
an optically thick medium, the polarization comes from 
upper layers. We then expect the field morphology to 
rather result from an upper surface layer than from the 
deeper embedded core. This is the first panel in Figure 
[7] which is labeled as envelope surface layer. The field is 
mostly uniform here, with a single orientation still reflect- 
ing the largest scale field structure. The field-to-gravity 
force ratio is relatively large (~ 0.7, Tabled]) and rather 
uniform with only a minor difference between core and 
outside-core region (Figure S). The correlation C, be- 
tween field and intensity gradient orientations, is poorest 
outside the core (C — 0.59). This all is suggestive of a 
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magnetic field that is still retaining its tension force in 
this surface layer. 

The JCMT 850 /jm observation - with a higher res- 
olution and reaching deeper inside the core - shows 
still some uniform patches in field morphology, but also 
clearly displays changing field orientations in the core. In 
this next deeper layer (labeled as molecular cloud in Fig- 
ure [7]) , the field ratio shows a clear difference between 
core (~ 0.4) and outside-core region (> 1, Figure Hf). 
This might indicate that gravity is overcoming the field 
tension in the core area, whereas the field remains sig- 
nificant at larger distances, possibly even holding back 
material from further accreting onto the core. Thus, the 
CSO and JCMT observations, despite covering similar 
areas, seem to sample different depths where the mag- 
netic field seems to play different roles. 

The SMA-SubC 870 /im observation is akin to zooming 
in on the JCMT observation. The force ratio Y>b and the 
correlation C are comparable to the JCMT values. This 
observation seems to further sharpen the JCMT core, 
removing some of the extended emission (third panel la- 
beled as core envelope in Figure [7|) . In addition to the 
JCMT observation (which already shows a core being 
dominated by gravity), the North- South symmetry in 
the SMA-SubC field morphology provides a hint for yet 
another role of the magnetic field: material can be chan- 
neled along the field lines onto an east- west axis. Finally, 
the SMA-Ext data with sub-arcsecond resolution clearly 
resolve the ~ 0.3pc structure into 4 cores (fourth panel 
labeled as core/disk in Figure [7]). Distinct morphologies 
are found for individual cores, with all cores consistently 
showing the tightest correlations (~ 0.9) and clearly low 
force ratios (~ 0.4). At this scale, an hourglass-like mor- 
phology is becoming visible. 

In summary, we thus propose a schematic scenario in 
Figure for the change in the B field, from the largest 
scale envelope surface layer down to the scale where 
a collapsing core is being revealed. The scenario de- 
scribed above emerges from visual impressions from the 
field morphologies, combined with a quantitative analy- 
sis from correlations and force ratios. 

Finally, it is interesting to remark that generally the 
correlation coefficients C are found to be larger in the 
cores than outside of the cores in all 4 scales and/or 
layers. The mean <ps (the difference between </>vi and 
0b) and its standard deviation in the cores tend to be 
smaller than the ones outside the cores. Both the anal- 
yses of C and (p$ suggest that 0b has a tighter corre- 
lation with 0vi in the cores. Similarly, the force ratio 
of the B field strength and the gravity (Sb) in these 4 
data sets also reveal a common trend: Sb tends to be 
smaller in the cores than outside the cores, suggesting 
that gravity is generally playing a more important role 
within the core s. As already quan tified and observed for 
other sources (IKoch et al.N2012bn . we can thus consis- 
tently confirm that the magnetic field significance sys- 
tematically decreases toward the cores for various phys- 
ical scales. 

5. CONCLUSION AND SUMMARY 

We present new high-resolution SMA dust polarization 
observations toward the molecular cloud W51 North. Be- 
yond the polarization imaging we also apply new analysis 
techniques providing insight into the role and importance 



of the magnetic field. Furthermore, we make the first at- 
tempt to present a magnetic field picture from the largest 
scales (3 pc) - from previously published CSO and JCMT 
data - to the SMA detected collapsing cores (~ 60 mpc). 
The key points are summarized in the following. 

1. Our new SMA high-resolution observations in the 
subcompact (SMA-SubC, angular resolution 8 syn 
of 4") and extended (SMA-Ext, sy „ of 0'.'7) array 
configurations clearly detect and resolve the dense 
structures in W51 North. With the observed ther- 
mal dust Stokes / continuum emission at a wave- 
length of 870 /urn, four dense cores, SMA1 to SMA4, 
are resolved in an east-west dust ridge. 

2. The polarized thermal continuum at 870 /im is de- 
tected and resolved with an angular resolution up 
to 0'.'7 for the first time with the SMA. The polar- 
ization percentage P ranges from 1 % to 4 %. P 
is found to decrease with higher intensity for both 
the SMA-SubC and SMA-Ext data. This trend is 
also largely identical to the larger scale CSO and 
JCMT polarization data. 

3. The observed magnetic field configuration on the 
~ 0.3 pc scale with the SMA-SubC shows some 
mirror-symmetry features from North to South. 
Here, the field lines might be channeling and align- 
ing material from both North and South. At a 
~ 60 mpc scale resolution observed with the SMA- 
Ext, dense cores are, indeed, detected along the 
east-west mirror symmetry axis. Smooth changes 
in field morphologies are found within individual 
cores, whereas large overall changes in position an- 
gles can occur in between the cores. 

4. A correlation between dust intensity gradient po- 
sition angles (0vi) and magnetic field position an- 
gles (0b) is found in the CSO, JCMT and both 
SMA data sets. This correlation is further analyzed 
quantitatively A systematically tighter correlation 
between 0vi and 0b is found in the cores, whereas 
at larger distances the correlation decreases. 

5. Magnetic field to gravity force ratio (Sb) maps are 
derived using the newly d eveloped polarizati on - 
intensity gradient method ([Koch et al.ll2012al ). In 
this method, Sb is quantified by measuring angles 
that reflect the resulting forces shaping the dust 
Stokes I and magnetic field morphologies. We find 
that the force ratios tend to be small (Sb ^ 0.5) in 
the cores in all 4 data sets. In regions outside of the 
cores, the ratios increase or the field is even domi- 
nating gravity (Sb > 1). This possibly provides a 
physical explanation of the tightening correlation 
between 0vi and 0b in the cores: the more the B 
field lines are dragged and aligned by gravity, the 
tighter the correlation is. Intuitively, in such a sce- 
nario one expects gravity to play an increasingly 
dominating role. The analysis via the force ratio 
Sb seems to consistently confirm this expectation. 
We stress that correlation and Sb are independent 
approaches. They, nevertheless, seem to reveal a 
matching picture. 
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Interpreting the four polarization observations 
covering different physical scales, we propose a 
schematic scenario for the magnetic field in W51 
North (Figure [7|). In the largest envelope sur- 
face layer (CSO), the field has kept much of its 
tension force. In the deeper-layer JCMT observa- 
tion, gravity has overcome the magnetic field in the 
core whereas in the outside the field tension is sig- 
nificant, capable of holding back material against 
gravity. Both higher-resolution SMA observations 
show gravity dominating the field tension. In SMA- 
SubC, the field is possibly channeling material from 



North and South onto an east-west axis. In SMA- 
Ext, four higher-resolution cores are aligned in an 
east- west axis, with a field morphology pointing 
toward hourglass-like structures as expected in a 
core/disk system. 
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